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The coherent flow of a Bose- Einstein condensate through a quantum dot in a magnetic waveguide is 
studied. By the numerical integration of the time-dependent Gross-Pitaevskii equation in presence 
of a source term, we simulate the propagation process of the condensate through a double barrier 
potential in the waveguide. We find that resonant transport is suppressed in interaction-induced 
regimes of bistability, where multiple scattering states exist at the same chemical potential and the 
same incident current. We demonstrate, however, that a temporal control of the external potential 
can be used to circumvent this limitation and to obtain enhanced transmission near the resonance 
on experimentally realistic time scales. 



PACS numbers: 03.75.Dg, 03.75.Kk, 42.65.Pc 

The rapid progress in the fabrication and manipula- 
tion of ultracold Bose-Einstein condensates has lead to a 
number of fascinating experiments probing complex con- 
densed matter phenomena in perfectly controllable en- 
vironments, such as the creation of vortex lattices [1] 
and the quantum phase transition from a superfiuid to 
a Mott insulator state in optical lattices [2]. With the 
development of "atom chips" [3-5] , new perspectives are 
opened also towards mesoscopic physics. The possibil- 
ity to generate atomic waveguides of arbitrary complex- 
ity above microfabricated surfaces does not only permit 
highly accurate matter- wave interference experiments [6] , 
but would also allow to study the interplay between in- 
teraction and transport with an unprecedented degree of 
control of the involved parameters. The connection to 
electronic mesoscopic physics was appreciated by Thy- 
wissen et al. [7] who proposed a generalization of Lan- 
dauer's theory of conductance [8] to the transport of non- 
interacting atoms through point contacts. Related theo- 
retical studies were focused on the adiabatic propagation 
of a Bose-Einstein condensate in presence of obstacles 
[9-12], the dynamics of soliton-like structures in waveg- 
uides (e.g. [13]), and the influence of optical lattices on 
transport (e.g. [14]), to mention just a few examples. 

Particularly interesting in this context is the propaga- 
tion of a Bose-Einstein condensate through a double bar- 
rier potential, which can be seen as a Fabry-Perot inter- 
ferometer for matter waves. In the context of atom chips, 
such a bosonic quantum dot could be realized by suit- 
able geometries of microfabricated wires. An alternative 
implementation based on optical lattices was suggested 
by Carusotto and La Rocca [15,16] who pointed out 
that the interaction-induced nonlincarity in the mean- 
field dynamics would lead to a bistability behaviour of 
the transmitted flux in the vicinity of resonances. This 
phenomenon is well known from nonlinear optics [17] and 
arises also in electronic transport through quantum wells 
(e.g. [18-20]) due to the Coulomb interaction in the well. 

In this Letter, we investigate to which extent resonant 
transport through such a double barrier potential can 



be achieved for an interacting condensate in a realistic 
propagation process, where the magnetic guide is grad- 
ually filled with matter wave. To simulate such a pro- 
cess, we numerically integrate the time-dependent Gross- 
Pitaevskii equation in presence of a source term that 
models the coupling to a reservoir of Bose-Einstein con- 
densed atoms. We shall point out that resonant scatter- 
ing states, which exist in principle for arbitrarily strong 
interactions, cannot be populated in the above way if the 
nonlincarity induces a bistability regime near the reso- 
nance. Finally, we suggest an adiabatic control scheme 
that permits to circumvent this limitation on experimen- 
tally feasible time scales. 

We consider a coherent beam of Bose-Einstein con- 
densed atoms propagating through a double barrier po- 
tential in a magnetic waveguide. In presence of a strong 
cylindrical confinement with trapping frequency luj_, the 
mean-field dynamics of the condensate is described by 
the effective one-dimensional Gross-Pitaevskii equation 
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with g — 2a s hto± [21], where m is the mass and a s the 
s-wave scattering length of the atoms. For the sake of 
defmitness, the double barrier potential is given by 



V(x) = V 
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(see Fig. 1). Our numerical calculations were performed 
for 87 Rb atoms (a s — 5.77nm) with uj± — 2ir x 10 3 s _1 , 
a± = y/h/moj± ~ 0.34/^m, Vq — hu>±, and L = 10a = 
5^m ~ 14.7o_l [22]. This yields g ~ 0.0347iu;xax- 

Let us first discuss resonances in terms of stationary 
scattering states of the condensate. The latter are given 
by stationary solutions ip(x,t) — ip(x) exp(— ifit/h) of 
Eq. (1) satisfying outgoing boundary conditions of the 
form ip(x) = Ae lkx with k > for x — ► oo. To calcu- 
late them, we insert the ansatz ip(x) — A(x) exp [i<j){x)\ 
(with real A and <f>) into the stationary Gross-Pitaevskii 
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FIG. 1. External longitudinal potential V in units of hui±. 
The dotted line shows the longitudinal atom density (in units 
of aj 1 ) of the scattering state associated with the 5th excited 
resonance, calculated at fj, — \.YlThus\_ and j t = 1.6u±. 



equation, and separate the latter into real and imagi- 
nary parts. This yields the condition that the current 
j(x) — (h/m)A 2 (x)<p' (x) = jt is independent of x, and 
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as equation for the amplitude A(x). The latter can be nu- 
merically integrated from the "downstream" (x — > oo) to 
the "upstream" (x — > —oo) region by means of a Runge- 
Kutta solver. As "asymptotic condition" at x — ► oo, we 
choose A' = Q and A satisfying 
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for a given jt, with n = A 2 the longitudinal density of 
the condensate. As was pointed out in Ref. [9], Eq. (4) 
exhibits a low-density (supersonic) and a high-density 
(subsonic) solution, where the transport is respectively 
dominated by the kinetic energy and by the mutual in- 
teraction of the atoms. Since in realistic propagation 
processes the waveguide is initially empty in the down- 
stream region, we choose the low-density solution for the 
asymptotic value of A. 

A measure of the proximity of the scattering state to 
a resonant state is provided by the drag 
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that the condensate exerts onto the obstacle [11]. Far 
from any resonance, the amount of reflection from the 
potential is rather large, leading to a finite drag due 
to the associated momentum transfer, while a vanishing 
drag is expected near a resonance where the condensate 
is perfectly transmitted through the quantum dot. In 
Fig. 2, the drag is plotted as a function of \i and j t in 
the vicinity of the 5th and 6th excited resonance, which 
have five and six nodes within the well, respectively (see 
Fig. 1). While the chemical potential of the resonant 
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FIG. 2. Drag exerted by the condensate onto the obsta- 
cle, plotted as a function of the chemical potential and the 
total current jt (in units of uj±) for g = (left panel) and 
g = 0.034:huj±a± (right panel). Light gray areas correspond 
to a low and dark gray areas to a high drag. In white areas, 
the integration of Eq. (3) leads to a diverging A(x). The lo- 
cation of the 5th and 6th excited resonance, where the drag 
vanishes, is marked by dashed lines. 



state is independent of j t in the linear case, it increases 
with jt in presence of a repulsive atom-atom interaction. 
Note that for large currents scattering states with non- 
diverging amplitude A(x) exist only in the immediated 
vicinity of resonances. 

Can resonant scattering states be populated in a realis- 
tic experiment where the condensate is initially confined 
in a microtrap and then released to propagate through 
the waveguide? To address this question, we numerically 
integrate the time-dependent Gross-Pitaevskii equation 
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in presence of an inhomogeneous source term emitting 
coherent matter waves with chemical potential \i at po- 
sition xq (we used xq = — Wa± in our calculation). The 
wavefunction ip(x,t) is expanded on a lattice (within 
—20 < x/aj_ < 20) and propagated in real time domain. 
To avoid artificial backscattering from the boundaries 
of the lattice, we impose absorbing boundary conditions 
which are particularly suited for transport problems [23] 
and can be generalized to account for weak or moderate 
nonlinearities [24,25]. 

Stationary scattering states can now be calculated by 
propagating ip(x,t) in presence of an adiabatic increase of 
the source amplitude So from up to a given maximum 
value, with the initial condition ip = 0. This approach 
simulates a realistic propagation process where a coher- 
ent beam of Bose-Einstein condensed atoms with chemi- 
cal potential y, is injected into the guide from a reservoir. 
Furthermore, it provides a straightforward access to the 
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FIG. 3. Transmission spectrum obtained from the 
time-dependent propagation approach for g — 0.034hiLi±a± 
and fixed incident current ji — 1.6ui± (solid line) compared, 
in the upper panel, to the interaction-free spectrum (dotted 
line) . The crosses show the results of a full three-dimensional 
calculation [24] . The dashed line in the lower panel shows the 
two other branches of the 5th resonance peak, which are not 
populated by the time-dependent propagation process. 

transmission coefficient T that is associated with a given 
scattering state: T can be defined by the ratio of the cur- 
rent j t in presence of the double barrier potential (i.e., 
the transmitted current) to the current ji obtained in ab- 
sence of the potential (the incident current). The latter is 
analytically evaluated as ji = h\So\ 2 / (mfc ) with k being 
self-consistently defined by k = ^/2m(tt — g\So\ 2 /k^j/fi. 

Fig. 3 shows the transmission coefficient as a func- 
tion of the chemical potential around the 5th excited 
resonance. For each value of /x, the wavefunction was 
propagated according to Eq. (6) in presence of an adi- 
abatic increase of the source amplitude S a up to the 
maximum value that corresponds to the incident current 
ji = 1.6lu±. The transmitted current is directly eval- 
uated from the stationary scattering state obtained at 
the end of the propagation procedure. While the typ- 
ical sequence of Breit-Wigner resonances is obtained in 
the linear case (or in the limit of very low incident cur- 
rents), the profiles become asymmetric for g > with 
perfect transmission being suppressed for narrow reso- 
nances. These results are essentially reproduced by a full 
three-dimensional mean-field calculation, which will be 
described elsewhere [24]. 

The step-like structures in the transmission spectrum 
indicate a bistability phenomenon, as known from similar 
processes in nonlinear optics [17] and in electronic trans- 
port through quantum wells (e.g. [18-20]). Additional 
branches of the resonance peaks are indeed identified by 



the integration method based on Eq. (3) which allows to 
calculate stationary scattering states for given j t and /i. 
The incident current of the scattering state is approxi- 
mately determined according to Ref. [12] via 



^ ("av + V n max™min) (7) 

with n av = 5(n max + n m i n ), where n max and n m i n denote 
the maxima and minima, respectively, of the longitudinal 
upstream density. The expression (7) assumes a cosine- 
like oscillation of the upstream density, and is valid for 
small g(n max - n m in)/n a v 

Finding the value of jt that results from a given ji is 
now an optimization problem that can be solved straight- 
forwardly. The result is shown in the lower panel of 
Fig. 3 for chemical potentials around the 5th resonant 
state. In addition to the spectrum obtained by the 
integration Eq. (6), two further solutions appear for 
1.02 < /j,/(hu)±) < 1.13 which join together to form a 
resonance peak that is asymmetrically distorted towards 
higher \i. The existence of such a multivalued spectrum, 
which is reminiscent of nonlinear oscillators, was in this 
context pointed out by Carusotto and La Rocca [15]. 
Since the additional branches of the resonance peak are 
apparently not populated by the time-dependent inte- 
gration approach, we expect that resonant transport will 
generally be suppressed in a realistic propagation pro- 
cess. Loosely speaking, the atoms "block" each other 
when going through the double barrier potential [26] . 

To enhance the transmission of matter waves near a 
narrow resonance, the external potential needs to be var- 
ied during the propagation process. Specifically, this can 
be achieved e.g. by illuminating the scattering region 
with a red-detuned laser pulse. We can describe such 
a process by a temporal modulation of V according to 

V(x)^V(x,t) = V(x)-V (t) (8) 

where Vo(t) > is determined by the detuning and the 
time-dependent intensity of the laser. For an adiabatic 
modulation of V, the wavefunction ip{x,t) will, at each 
time t, remain close to the instantaneous scattering state 
associated with the external potential (8) — or, equiva- 
lently formulated, to the scattering state for V = V(x) 
at the shifted chemical potential /i + Vo (t) . Outside the 
bistability regime, e.g. for /j, + Vo(i) < 1.02huj± in case of 
the 5th resonance, this scattering state would be uniquely 
given by the one that is also obtained by the direct propa- 
gation process. However, as soon as the effective chemical 
potential /x + Vo(i) is raised above 1.02 Huj±, the wave- 
function will continuously follow the upper branch of the 
resonance and evolve into a near-resonant scattering state 
with high transmission. 

Indeed, we can use our numerical setup to simulate 
such a process. Fig. 4 shows the transmission coefficient 
as a function of the propagation time, where the effective 
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FIG. 4. Time evolution of the transmission coefficient T 
(upper panel) during the ramping process of the external po- 
tential Vo (lower panel) which shifts the effective chemical 
potential from /i = 0.985 to 1.1257ia>x- As shown in the in- 
sets (with scaling as in Fig. I), the wavefunction adiabatically 
evolves into a nearly resonant state with transmission close to 
unity, and decays from there to the low-transmission scatter- 
ing state within a time scale of the order of r ~ lOOwJ 1 . 



chemical potential was shifted from /i = 0.985 TujJ±_ to 
/i + V maiK = 1.125 fiui± by means of a Gaussian ramping 
process taking place within < uj±t < 360 (see the lower 
panel of Fig. 4). We see that the transmission approaches 
unity at the end of the ramping process, which clearly 
indicates that the scattering wavefunction evolves along 
the upper branch of the resonance peak. This is indeed 
confirmed by the associated density shown in the insets 
(to be compared with Fig. 1). 

As was also pointed out in the context of electronic 
transport through quantum wells [20] , the resonant scat- 
tering state is dynamically unstable in presence of in- 
teractions. This instability is indeed encountered in 
our system: Continuing the numerical propagation be- 
yond uj±t = 360 (at fixed Vo) results in a decay of the 
wavefunction towards a low-transmission scattering state 
within a time scale of the order of r ~ lOOwJ 1 ~ 16ms. 
This lifetime should be long enough, however, to trans- 
port a large fraction of the condensate through the double 
barrier, as well as to manipulate the resonant scattering 
state: by closing, for instance, the potential well during 
that time scale (e.g. with blue-detuned lasers that en- 
hance the potential outside the barriers), one would cre- 
ate a trap in which an interacting mean-field state with 
an unusually high excitation (with five nodes in case of 
the 5th excited resonance) would be obtained. 

In conclusion, we have studied resonant transport of in- 
teracting Bose-Einstein condensates through a symmet- 
ric double barrier potential in a magnetic waveguide. The 



nonlinearity induced by the interaction leads to a distor- 
tion of the resonance peak, where the associated scatter- 
ing state cannot be populated by directly sending coher- 
ent matter wave onto the initially empty waveguide. To 
obtain nevertheless a finite amount of transmission on 
intermediate time scales, the external potential needs to 
be adiabatically varied during the propagation process. 
The lifetime of the resonant scattering state obtained in 
this way is predicted to be of the order of r <~ 10ms for 
our particular setup, which should be long enough to al- 
low for further experimental manipulations. We expect 
that the basic principles of the scenario encountered for 
our double barrier potential apply also to more complex 
quantum dot geometries such as sequences of more than 
two barriers along the guide. This indicates that the de- 
sign of suitable control schemes will be a relevant issue for 
the mesoscopic transport of Bose-Einstein condensates. 

It is a pleasure to thank Nicolas Pavloff, Peter 
Schmclchcr, Joachim Brand, Jozsef Fortagh, and Wil- 
helm Prcttl for fruitful and inspiring discussions. 



[1] J. R. Abo-Shaeer, C. Raman, J. M. Vogels, and W. Ket- 

terle, Science 292, 476 (2001). 
[2] M. Grciner et al, Nature 415, 39 (2002). 
[3] R. Folman et al, Phys. Rev. Lett. 84, 4749 (2000). 
[4] H. Ott et al, Phys. Rev. Lett. 87, 230401 (2001). 
[5] W. Hansel, P. Hommelhoff, T. W. Hansen, and J. Reichel, 

Nature 413, 498 (2001). 
[6] E. Andersson et al, Phys. Rev. Lett. 88, 100401 (2002). 
[7] J. H. Thywissen, R. M. Westervelt, and M. Prentiss, 

Phys. Rev. Lett. 83, 3762 (1999). 
[8] D. K. Ferry and S. M. Goodnick, Transport in Nanostruc- 

tures (Cambridge University Press, Cambridge, 1997). 
[9] P. Leboeuf and N. Pavloff, Phys. Rev. A 64, 033602 

(2001). 

[10] M. Jaaskelainen and S. Stenholm, Phys. Rev. A 66, 

023608 (2002). 
[11] N. Pavloff, Phys. Rev. A 66, 013610 (2002). 
[12] P. Leboeuf, N. Pavloff, and S. Sinha, Phys. Rev. A 68, 

063608 (2003). 

[13] S. Komineas and N. Papanicolaou, Phys. Rev. Lett. 89, 
070402 (2002). 

[14] K. M. Hilligs0e, M. K. Oberthaler, and K.-P. Marzlin, 

Phys. Rev. A 66, 063605 (2002). 
[15] I. Carusotto and G. C. La Rocca, Phys. Rev. Lett. 84, 

399 (1999). 

[16] I. Carusotto, Phys. Rev. A 63, 023610 (2001). 
[17] R. W. Boyd, Nonlinear Optics (Academic Press, London, 
1992). 

[18] V. J. Goldman, D. C. Tsui, and J. E. Cunningham, Phys. 

Rev. Lett. 58, 1256 (1987). 
[19] C. Presilla, G. Jona-Lasinio, and F. Capasso, Phys. Rev. 

B 43, 5200 (1991). 
[20] M. Y. Azbel', Phys. Rev. B 59, 8049 (1999). 
[21] M. Olshanii, Phys. Rev. Lett. 81, 938 (1998). 
[22] Such a small barrier width would be at the limit of real- 



4 



izability with present-day atom-chip technology. We re- 
mark, however, that the phenomena we discuss here do 
not sensitively depend on the chosen parameters nor on 
the specific shape of the potential. 

[23] T. Shibata, Phys. Rev. B 43, 6760 (1991). 

[24] T. Paul, P. Schlagheck, and K. Richter, in preparation. 

[25] We reproduced our results by an alternative approach 
using complex absorbing potentials at the boundaries. 

[26] This phenomenon must not be confused with the "atom 
blockade" discussed by Carusotto [16], which arises in the 
limit of few atoms and very strong interactions. 



